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A  Maximum-Likelihood  Based  frequency 
Synchronizer  for  Dual-h,  Pull  Response  4-ary  CPM 

Ray  H.  Pettit  and  Bruce  E.  Wahlen 


Abstract 

Maximum-likelihood  techniques  are  useful  in  finding  synchronizer  structures  for  various  cases.  Synchronizers  for  frequency,  phase,  and 
timing  have  been  found  for  various  bandpass  signaling  techniques  such  as  PSK,  DPSK,  QAM,  MSK,  and  CPM.  These  include  data-aided, 
decision-directed,  and  clock-aided  cases.  For  CPM,  however,  apparently  only  the  single  modulation  index  case  has  ML-based  synchronizers. 

This  report  describes  a  new  non-data-aided,  non-decision-directed  ML-based  frequency  synchronizer  (with  no  phase  or  timing  information), 
derived  for  a  full- response,  dual-h  (two  modulation  indexes),  4- ary  CPM  signaling  scheme.  The  derived  structure  will  be  incorporated  into 
future  simulations  to  compare  performance  among  several  possible  frequency  synchronizers. 
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L  INTRODUCTION 

Successful  operation  of  digital  communication  systems  requires  that  receivers  achieve  synchronization.  Carrier  fre¬ 
quency  and  phase,  waveform  time-of-arrival,  frame  synchronization  for  TDM/TDMA,  and  timing  for  a  frequency-hopping 
pattern  or  for  a  direct  sequence  spreading  code  for  spread-spectrum  are  examples  of  typical  unknowns  that  must  be 
estimated  for  good  performance. 

There  are  many  successful  synchronization  schemes.  Mengali  and  D’Andrea  [1]  present  theoretical  and  practical  details 
for  various  synchronizers  for  frequency,  phase,  and  symbol  timing.  Their  text  deals  with  baseband  and  various  bandpass 
signaling  techniques,  including  PSK,  DPSK,  QAM,  MSK,  and  CPM  (single  modulation  index  case).  This  p)aper  concerns 
an  important  special  version  of  CPM  (Continuous  Phase  Modulation),  namely  dual-h,  full-response  with  a  rectangular 
frequency  response  function,  and  4-ary  signaling  alphabet  [2].  In  addition  to  frequency,  phase,  and  symbol  timing,  this 
case  also  requires  the  so-called  “super-baud”  synchronization  inherent  with  CPM  using  several  modulation  indexes.  The 
paper  presents  a  structure  for  frequency  acquisition  and  tracking  for  this  case,  based  upon  an  optimum  procedure  known 
as  “maximum  likelihood”  [3].  This  is  a  standard  criterion  for  optimality  (maximization  of  a  probability  density  function 
for  received  samples  conditioned  upon  a  set  of  parameters  to  be  estimated).  Future  papers  will  address  other  levels  of 
synchronization  for  this  case. 

Synchronizers  often  develop  from  good  intuitive  reasoning  rather  than  from  estimation  theory.  (In  many  cases,  these 
ad  hoc  synchronizers  fit  within  an  estimation  theory  fi-amework.)  Others  have  been  invented  by  adopting  the  techniques 
of  estimation  theory  from  the  start.  Both  procedures  are  very  useful.  In  this  paper,  the  latter  approach  is  followed. 
Even  though  the  initial  modeling  and  derivation  efforts  are  for  this  “optimum”  maximum-likelihood  criterion,  it  should 
be  emphasized  that  the  inevitable  necessary  approximations  for  analytical  completion  and  receiver  implementation  lead 
to  sub-optimal  final  frequency  recovery  schemes.  Final  questions  involving  performance  of  this  sub-optimal  synchronizer 
will  require  simulation  methods  because  of  inherent  analytical  complexities.  A  direct  measure  of  the  “goodness”  of  an 
estimator  is  its  variance,  usually  compared  to  either  the  Cramer-Rao  or  the  Modified  Cramer-Rao  lower  bounds  (see  [1], 
for  example).  Ultimately,  desirable  performance  would  be  in  terms  of  effects  on  ability  to  correctly  receive  data  when 
using  the  particular  synchronization  method.  A  future  paper  will  present  results  from  appropriate  simulation  studies. 

II.  General  Structure  for  the  Frequency  Synchronizer 

The  task  at  hand  is  to  find  a  frequency  recovery  scheme  (acquisition  and  tracking)  that  can  operate  without  knowledge 
of  the  carrier’s  phase,  symbol  time-of-arrival,  or  data.  That  is,  the  synchronizer  is  to  be  non-data  aided  and  non-decision 
directed.  The  notation  and  modeling  closely  follow  Mengali  and  D’Andrea  [1]  as  used  for  the  single  modulation  index 
CPM  considered  in  their  text.  Complex  envelope  notation  is  used  for  convenience.  The  i>aper’s  primary  notation  is 
given  below. 

We  assume  that  the  received  waveform,  r(t),  consists  of  the  CPM  signal  and  an  additive  white  Gaussian  noise 
component  of  spectral  density,  ATo, 

r(f)  =  s(f)+^t)  (1) 

with 

s(f)  =  ,  (2) 

First  author:  Professor,  Electrical  and  Computer  Engineering  Department,  California  State  University,  Northridge,  CA  913S0-8346. 
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In  (2),  V,  and  r  represent  the  unknown  parameters:  frequency  offset  from  the  nominal  carrier  value,  carrier  phase,  and 
time-of-arrival,  respectively.  Es  is  the  “energy”  of  the  signal  waveform  over  the  symbol  interval,  T,  The  data  sequence 
is  a  =  (•  •  • ,  a~i ,  ao,  ai ,  •  •  •),  where  ai  €  {±1,  ±3}  for  the  4-ary  case  of  interest. 

The  phase  a)  is  given  by 

cc)  =  27r  ^  aihiq(t  —  iT)  ,  (3) 

i 

where  hi  is  the  “modulation  index”  for  the  interval  For  our  dual-h  case,  two  modulation  indexes  are  used,  and 
hi ,  with  fiQ  associated  with  the  even  numbered  intervals  and  hi  with  the  odd  ones.  Thus, 


i){t,cx)  -  ^i(t,Q:e) (4) 
=  27r<ho^2 ^ ct2i+iq{t  -  [2i  -f-  1]T) 

It  t 

In  (3),  q{t)  is  the  phase  response  function,  which  is  shown  below  for  the  full-response  rectangular-frequency-pulse  case 
(so-called  1  REC): 

{0,  ift  <0 

ifo<i<r  (5) 

I  ift>T. 

In  keeping  with  modern  practice,  processing  of  discrete-time  samples  will  be  performed  by  the  synchronizer,  a  so-called 
“digital  synchronizer,”  instead  of  an  “analog  synchronizer”  that  processes  directly  the  continuous-time  waveforms.  See 
figure  1.  The  samples  x  =  are  taken  from  the  output  of  an  “anti-aliasing’  filter,  AAF,  which  is  assumed  to 

have  sufficiently  wide  passband  to  pass  through  the  signal,  ${t)^  without  significant  distortion.  The  sampling  period,  Tg, 
is  sufficiently  short  to  allow  the  samples  to  retain  the  information  content  of  the  continuous-time  waveforms. 


Fig.  1.  Sampling  output  of  anti-aliasing  filter. 


For  some  additional  notation,  we  let  Lo  be  the  number  of  information  symbols  to  be  “observed”  for  purposes  of 
collecting  the  samples,  and  let  N  be  the  number  of  samples  taken  for  each  symbol.  (Mengali  and  D’Andrea  suggest 
jV  =  4  as  an  appropriate  value  for  cases  studied  by  them.)  Thus, 

To  =  “observation  time”  (6) 

=  LoT 
=  NLoTs 


and 

X  =  {x(0),  x(7;), .  •  • ,  x{[NLo  -  1]T,)}  .  (7) 

Since  the  task  is  to  derive  a  maximum-likelihood  based  estimator  of  the  unknown  frequency  offeet,  v,  which  does 
not  require  knowledge  of  {0,r,  a},  we  require  the  “Likelihood  Function,”  A(x|S^,  which  will  be  obtained  from  the 
“Conditional  Likelihood  Function,”  A(a:|5,t;,  6,t),  by  averaging  over  {S,0,t}.  (See,  for  example.  Van  Trees  [3]  for  an 
excellent  treatment  of  detection  and  estimation  theory.) 

Before  continuing,  some  clarification  of  notation  is  in  order.  For  example, 

V  =  actual,  but  unknown,  frequency  offeet 
S’  =  a  possible  value  of  v  (a  “realization”) 

S  =  an  estimate  of  v  . 


The  conditional  likelihood  function  is  given  by 

T, 


A(a:|5,S,0,r)  =exp^  —  Re 


t 


NLo-t 

xikTs)r{kT,) 

L  fc=o 


2No 


E  iw)P 


A:=0 


(8) 
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5(i)  =  ,  (9) 

Since  the  last  sum  in  (8)  is  independent  of  |a,u,  0,t|-,  it  is  sufRcient  to  consider 

{j,  rivi-o-i  1 1 

^Re  ^  x{kTs)T{kT,)  |  (10) 

{rp  iOJP  ^NLq  —  I  ^ 

e"-’®  ^  x{kTs)e~^^'^'’'‘'^‘ 

L  A;=0 

By  definition,  the  “maximum  likelihood  estimate  of  v"'  is  the  value  oft;  that  maximizes  A(a;|i)), 

To  proceed,  we  define 

fc=0 

(the  magnitude  of  the  complex  X  is  |J\r|  and  its  angle  is  <j>x]  i.e.,  X  ^\X\  ),  Noting  that  X  is  independent  of  0,  we 

have 

A'(x|a,?;,?,r)  =  expIciXlcos  (4>x  ~^)}  (12) 

with 

^  %  1^ 

T" 

First  we  average  over  0,  assuming  that  9  is  uniform  (0,27r): 


A'(x|q,z>,t)  =  -^  J  exp|c|X|cos  d?| 

=  /o(C'|X|). 


/o(*)  is  the  zeroth  order,  modified  Bessel  Function  of  the  First  Kind,  a  real  monotonically  increasing  function  of  its 
argument. 

It  is  not  analytically  feasible  to  obtain  the  exact  Vji^i  at  this  point,  so  we  proceed  based  on  the  special  case  of  low 
signal-to-noise  ratio.  This  allows  the  approximation, 

I^{C\X\)^l  +  ^\Xf  ,  (14) 

and  leads  ultimately  to  a  maximum-likelihood  based  estimator  of  v.  Maximizing  A'(a;|?5)  is  therefore  approximately 
equivalent  to  maximizing 


A"(x\^  =  E^~{\Xf‘}  , 


where  (15)  implies  averaging  \xf  over  {cx,T).  Appendix  A  presents  the  details  of  this  averaging,  showing  that 

NLo-lNLo-l  ^ 

A"(x|tf)=  £  £  x(feiT,)x*(fc27;)e--'N‘‘-**)^^*]/f[fci7;,A:2r,]  , 

ki=0  *2=0 


where  H  [kiTsyk2Ts]  is  defined  by 


-J  ^2T  oo  r 

H[k^T„k2Ts]  =  ^  n 

t=:-00  • 


-  TT  1  sin(8^ftoP  %Ts  -  T  -  2iT,  (kj  -  ki)T,]) 
^  Jo  i=-«  L ^  sin(27r/iop  [*27;  -  r  -  2iT,  (*2  -  *i)7;j) 

1  sin(87rAip  fer.  -  f  -  (2i  +  1)T,  (*2  -  ki)T,])  1 
4  sin(27rAip  [*2?;  -  r  -  (2i  +  1)T,  (*2  —  *1 )  Ji])  ^  ’ 


In  (17),  p(t,  At)  is  defined  by 


p{t,At)  =  q{t)-q(t- At)  . 
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By  changing  the  variable  of  integration  in  (17)  to  t  =  k2Ts  —  r  and  noting  that  the  integrand  is  periodic  with  period 
2T,  we  have 


H[k^Ts,k2Ts] 


H[{k2-ki)Ts]^H[{ki-k2)Ts] 

1  -fj  risin(87rftop[t-2iT,(A^-A:i)r,]) 

27^  Jo  [4  sin(27r/ioP  [t  -  2iT,  (A^  -  feijr,]) 


1  sin(87r/iip  \t  —  {2i  +  l)r,  (A^  —  A:i)Ti,]) 

4  sin(27r/iip  [t  —  (2i  4-  \)T,  (k^  —  A:i)r4]) 


dt  . 


(19) 


The  general  structure  for  the  frequency  synchronizer  can  now  be  developed,  paralleling  Mengali  and  D’Andrea’s 
development  for  the  single-h  case  [1],  although  significant  additional  work  must  be  carried  out  to  produce  a  synchronizer 
of  reasonable  implementation  complexity.  Section  III  will  describe  the  additional  work. 

We  seek  the  value  of  v  that  maximizes  A"(a;|v),  as  shown  in  (16).  We  take  the  derivative  and  set  the  result  to  zero, 
in  the  usual  way: 


iVLo-l  JVLo-l 

=  i27rT.  £  £  (20) 

fcl=0  ^2=0 

•  -  h)H  [(*2  -  ki)Z]] 

=  0. 


Because  A"(x|v)  is  real,  so  is  its  derivative,  and  thus  the  term  in  (20)  represented  by  the  double  summation  is  imaginary. 
Now  define 


y{kTs)  = 

h{kTs)  =  kH{kTs)  . 


We  seek  the  value  of  v  to  satisfy 


Let 


{NLo-l 

£  yikiT,) 

A:i=0 


NLo-l 

y{k2T,)h[{ki-k2)T,] 

.  fc2=0 


=  0 


z{kTs) 


NLo-l 

£  yik2Ts)h[ik-k2)T,] 
*2=0 


y{kTs)  *  h(kn)  . 


(21) 


(22) 


(23) 


(Since  h(kTs)  is  real,  z^kZ)  =  Zk,i°o~"  V^iW^k  -  A^)T.]). 

Thus,  h{kTs)  can  be  interpreted  as  a  digital  filter’s  impulse  response  and  z{kTs)  as  its  output  with  y{kTg)  being  the 
input. 


y(kZ) 

. 

Filter 

z(kZ) 

hikZ) 

- K 

Before  proceeding  to  the  final  general  structure,  we  recognize  that  h{kTg)  is  non-causal.  This  requires  that  a  suitable 
delay  of  D  sampling  intervals  be  incorporated  into  h{kTs)  to  have  physical  realizability.  Mengali  and  D’Andrea  [1]  for  a 
similar  situation  state  that  D  =  2N  is  probably  adequate,  but  further  study  of  this  point  will  be  needed  for  our  case  and 
will  be  i>art  of  the  associated  simulation  studies.  In  Section  III  we  will  provide  considerably  more  detail  about  ^(AjTfi), 
thus  providing  the  basis  for  the  final  simulation  steps. 
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We  show  the  delay,  D,  as  follows: 


wikT,)=i{k-D)T,] 


Now  we  re-write  (22)  as  follows: 


NLq~\+D 

Y,  {y  [(A  -  D)T,]  w\kTs)]  =  0  .  (24) 

fci=D 

Equation  (24)  provides  the  basis  for  the  general  structure.  The  basic  idea  is  to  exploit  the  sum  of  some  consecutive 
terms  in  the  above,  say,  N  terms,  as  an  error  signal  to  drive  — towards  zero. 

dv 

The  frequency  estimates  are  to  be  updated  according  to 

V  [(n  +l)T]=v  [nT]  +  'ye  [nT]  ,  (25) 

where  7  is  a  step-size  parameter  and  e  [nT]  is  the  error  signal, 

{n+l)N-\ 

e[nTl=  Y  lra{y{kTs)w\kT,)}  .  (26) 

k=nN 

The  result  is  the  structure  shown  in  Figure  2.  (In  (26),  n  is  the  symbol  index  and  k  is  the  sample  index). 

The  approach  described  above  is  a  standard  one  for  finding  synchronizer  structures  with  closed-loop  feedback  process¬ 
ing.  (See  [l].)  The  parameter,  7,  selected  based  on  simulation  studies,  should  not  be  so  large  as  to  cause  large  variations 
of  the  estimate  about  the  true  value  nor  so  small  as  to  require  excessive  time  to  converge  to  the  true  value. 


<[(n  +  V)T]=v(nT)  +  ye(nT) 


Fig.  2.  General  structure  for  frequency  synchronizer 

The  general  structure  shown  in  Figure  2  is  identical  to  Figure  4.10  of  [1]  for  the  single-h  CPM  case.  The  effects  of 
dual-h  CPM  are  incorporated  in  the  details  of  /i(A:T5),  which  seem  to  be  significantly  different  from  the  h{kTs)  of  single-h 
CPM  in  [1].  Section  III  will  now  derive  and  discuss  these  additional  details. 

III.  The  Feedforward  Filter,  h{kTs) 

The  function,  p{t,kTs)  (and  all  versions  of  it  shifted  by  integer  multiples  of  NTg)^  holds  the  key  to  understanding 
H  [(^2  ^i)Js]5  shown  in  (19),  and  to  implementing  h{kTs)  =  kH  Because  p(^,  kTs)  h2is  finite  width  (see  Appendix 

B  for  details  of  p{t,kTs)  as  a  function  oi  k  =  k2  —  ki),  there  are  only  a  finite  number  of  factors  in  the  product  of  (19) 
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that  differ  from  unity.  That  is,  if  both  p\t  —  2iT,  kTg]  and  p[t—  {2i  4-  1)T,  kTg]  are  zero  over  the  range  of  the  integral, 
(0, 2T)  (for  a  particular  value  of  k),  that  particular  ith  factor  equals  unity.  Therefore,  the  only  values  of  i  in  the  product 
of  (19)  that  are  significant  to  the  understanding  of  the  behavior  of  H{kT8)  are  those  for  which  the  corresponding  factors 
differ  from  unity. 

In  Appendix  C,  these  “significant”  values  for  i  are  shown  as  follows: 


k 

Set  of  Significant 
Values  of  i 
(from  Equation  (19)) 

0<k^  N 

0,-1 

k>N 

-N^k<0 

0,1 

A 

1 

In  the  above,  represents  the  next  integer  equal  to  or  exceeding 

Using  the  above,  N\x\v)  from  (16)  can  be  re-written  as  follows,  requiring  only  a  finite  number  of  factors  to  be 
considered: 

NLo-l  _ 

M'{x\v)=  2  (27) 

where 


and 


ki-N-1 

5 

k2 


r2T  L  J  int 


^  £  n  (•) 

*=1-1  „  ,  f2T  1 

k^^-N  -^0  r=o 

V+fei  _  .2T  0 

=^+i  Jo 


dt 


fc2  =  fcl  +  l 
NLo-1 


k2=N+ki+l  •'0  . _ rfc,-ti  1H 

*-  1  2iV  Ji 


{•)dt 


1  sin(87r/iop  [t  —  2iT,  {k^  —  A:i)7i]) 

4  sin(27rftoP  [*  —  2iT,  (Aa  —  *i)7»]) 

1  sin(87rfeip  [t  —  (2i  -f-  1)T,  (fca  —  fei)Ts]) 
4  sin(27r/iip  [i  —  (2i  +  1)7’,  (/ca  —  *i)?*]) 


(28) 


(29) 


(See  (19)). 

The  next  step  for  understanding  and  implementing  h{kTs)  is  the  evaluation  of  the  four  integrals  of  (28).  Many  tedious 
steps  are  required.  From  Appendix  D,  we  have  the  following  results: 

Case  1:  0  <  k2  <  ki N  —  I 


Fi  (A:i!7^,  k2Ts) 


r2Tl  \ir.t 


1  rZl  I.  J  int 

-  n  (•)- 

,  1 1  r»=i-*2]+ 

=  [cosTr/iocosvrfti  cos27r/iocos27rfti]  I  !*»< 

=  Fi[(*2-A:i)T,] 


(30) 
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Case  2:  ki-N<k2<ki  —  l 


^2(,kiTg,k2Ts)  = 


2T 

i=0 

(^2  ~  fel) 
N 


3nho(k2  -  fci)  .  Trho{k2  —  ki) 

cos  -  +COS  - 


,  Stt/ii  (k2  —  ki)  IT  hi  ik2  —  ki) 

+  COS - - - +COS - i 


+ 

4- 

+ 

+ 

F2 


Iff  hi 

477  \  [3(/i2  -  hj 


+ 


Bhi 


hi 


hf)  h^  —  Qh\ 

hi 


.  Z-Khi{k2  —  ki) 
- N - 


Qh^-hl  hi-  h\ 

ho  ,  ho 


[9hi-hl^  hi- hl\ 
ho  ,  Zho 


Sin 


sin 


3(ftf  -  hi)  hi-  9hl\ 
{k2  -  ki)T,] 


'Khi{k2  —  ki) 

N 

irho{k2  —  ki) 

N 

.  37r/io(&2  —  *i) 
- N - 


Case  3:  ki+1  <k2<N+ki 


(31) 


1  ° 

F3ikiT.,k2Ts)  =  — jf  JJ  (.) 


dt 


=  1  [i  (fe2-fcl) 

4  N 


Z-!rho(k2  —  fci)  Tr/ioCfo  —  ^i) 

oas - - - +COS - - - 


Snhiik^  —  ki)  Trhiiko  —  kA 

cos - - - +COS - - - 


+ 

+ 

+ 

+ 

=  F3 

Case  4:  N  +  ki  +  l<k2<  NLo  —  1 


j_  rr _ ^ 

477 


Zhi 


[hi  -  9hl  hi- hi 
ho  .  ho 


hi 


hl)^9hl-hl\ 

hi 


sin 


ZTrhi{k2  —  ki) 


N 


Ag  -  9hl  hi-  hl\ 


ho 


Z{hl  —  hi)  ^  9hl  —  hl\ 
{k2-ki)Ts] 


Zho 


.  irhi(k2  —  ki) 

N 

.  Tvho{k2  -  All) 

TV 

.  Z'Kho{k2  —  ki) 
TV 


(32) 


r2T 


F^ikiTsMTs)  =  ^  /  n  (33) 

ir=  [>=^-*.1  + 

L  2"  Ji„« 

-  ,  r  ~ fci  1 

=  [costt/iq  costt/ii  cos27r/io  cos27r/ii]  iirtt 

=  F4[{k2-’ki)Ts] 

As  discussed  in  Appendix  D,  analytical  complexities  led  to  the  need  for  the  approximations  used  to  obtain  Fy  [(/C2  — 
and  F4  [(A:2  -  A:i)Ts].  Exact  results  were  obtained  for  F2  [(A:2  -  A:i)T5]  and  F3  [{k2  -  ki)Ts]- 

It  appears  that  Fj  [(A:2  —  A)i)r5]  andF4  [(^2  —  make  only  small  contributions  to  G{ki,k2,^  in  (28)  and,  therefore, 

they  will  be  discarded  for  the  sequel.  This  allows  for  a  simple  implementation  of  h{kTs)  and  provides  for  a  reasonable 
delay,  D.  Note  that  F2  [(/:2  —  ^1)^5]  and  F3  [(/C2  —  fci)r5]  are  easily  calculated  for  fixed  values  of  ho  and  hi. 


8 


The  digital  filter,  h{kTs)^  in  Figure  2  can  now  be  more  specifically  given  by 

h{kTs)  =  kH{kTs) 

=  kF2{kTs)  ,  for  -  iV  <  A;  <  --1 
=  kFz{kTs)  ,  for  1  <  /:  <  AT 
=  0  ,  for  A;  =  0  and  [^1  >  AT . 

(See  Equations  (31),  (32)  for  F2{kTa),  Fs{kTs)^  respectively).  Note  that  H  [AjTs]  is  an  even  function  of  k  and,  therefore, 
^(A:T5)  is  odd. 

IV.  Summary  and  Conclusions 

This  paper  presents  the  derivation  of  a  frequency  synchronizer  for  dual-h,  full-response,  4-ary  CPM  (rectangular 
frequency  response  function).  The  synchronizer  is  based  upon  the  Maximum-Likelihood  criterion,  with  the  final  structure 
resulting  from  certain  approximations  made  to  overcome  various  mathematical  complexities  and  to  make  any  subsequent 
implementation  of  reasonable  complexity. 

It  is  not  clear  that  this  new  synchronizer  provides  significant,  if  any,  benefits  compared  with  others  in  use  or  proposed 
for  specific  applications.  Comparisons  based  upon  simulation  are  required  before  any  further  conclusions  can  be  drawn. 
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APPENDIX  A:  Derivation  of  A"(x|tJ) 

Since  A"(a;|i^  =  we  begin  with  (using  (11)) 

NLo-lNLo-l  _ 

ki  =0 

g+jVi  {*2^5  -'r,Q:c) 

(^1  (fc2  T»  -  TjSo) 


Because  the  data  symbols  are  independent,  the  averaging  over  cc  can  be  separated  into  the  Sg  and  5^  terms.  Consider 
the  first  term. 


•-r,ae)+J>i  -r,ae)  | 

(A2) 

=  F>-  «2i[9(A:iT,~7--2iT)-g(fc2T,-7-.2ir)])  | 

“  |gi27rAoa2i[9(fc2T,-r-2iT)"-9(/:ir,-r-2iT)]| 

i 

The  last  step  follows  because  the  oc2i  are  independent. 

By  defining 

p{t,At)  =  q{t)  -  q{t  -  At)  , 

(A3) 

the  above  product  can  be  rewritten  as 

|gi27r;voa2ip[A:2T;-T-2tT,(fc2“’^i)T>] 

i 

(A4) 

Now,  we  consider  each  factor  separately. 

^  1  gj  27r  /lo  a2i  P  [^2  Ts-t- 2iT,{k2  -  fci  )T^  ]  | 

=  i  |gi67r/iop[A:27;-7-2tT,(fc2-*i)T,] 

(A5) 

_|_  ^2TThop[k2Ts-7-2iT,{k2-ki)Ts] 

_j_  ^’-j2'Khop[k2 Ts  -7-  2iT,(k2  - ki  )T,] 

^  g“j67r;iop[fc2T,-7-2tT,(/!:2-fci)T^]|  ^ 
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which  is  of  the  form, 

-  -f  i  (cos 3/? 4-  cos /?) 

=  cos  /?  cos  2^ 

_  1  sin  2/3  cos  2/3 

2  sin/3 
_  1  sin  4/3 

4  sin  /3 

Thus, 

|gj27r/ioa2ip[fc2T,-T-2»T,(/:2““fci)r,]| 

_  1  sin  (Stt/iop  -  r  -  2zT,  {k2  -  ArQTs]) 

4  sin  {2'Khop  [k2Ts  -  r  -  2iT,  (^2  -  A:i)Ts]) 

Similarly, 

E'^  /gi27r/iiQ:2t+ip[fe2T,-T-(2i+l)T,(fc2-fci)T,]  1 

OC2i+l  \  / 

_  1  sin  (Stt/iip  [fe2T5  -  r  -  (22  +  1)T,  {k2  -  ki)Ts]) 

4  sin  (27rhip  [A^T^  -  r-(2i-h  1)T,  (A;2  -  A;i)T;]) 

Using  (A7)  and  (A8),  and  then  averaging  over  the  uniform  r,  (0,2T),  gives  (17). 

APPENDIX  B:  The  Function  p[t^kTs] 

The  function,  p  [i, kTg]  =  q{t)  —  q{t  —  AiT^),  is  a  function  of  k  and  is  shown  below  for  reference. 


9(0 


(A6) 


(A7) 


(A8) 


Table  of  At-T,kTs]  =  At-NT„kT,] 


At-NT^fi]=^  for  all  t 


0<k<N 


t  _  1 
2NT,  2 


k  t 
aN~7m. 


NT,  iN  +  k)T,  2NT,  (2N+k)T, 


2  2NT 


2NT,^^^^ 

1  1 

1  1 

1  1 

NT, 

2NT,  (N+k)T, 

(.2N  +  k)/,‘ 

(N  +  k)T,  NT, 

(2N  +  k)T,  2NT, 

N<k<0 

“1/  , 

2N  2  2NtW 

'/  2NT.  ^ 

{N  +  k)T,  (2N  +  k)T,  NT,  2NT, 


2N  2  2NT,  \ i 


1/  2NT, 


(note:  not  drawn  to  scale) 
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APPENDIX  C:  Determination  of  Significant  Values  of  i  (from  Equation  (19)) 

Case  1:  0  <  k<N 
For  this  set  of  k  values, 


p  [t  -  2iT,  kTs]  =  0  ,  for  2iNTs  >  2NTs 
{0  <  t<  2NTs)  ,  (t  >  1) 

and 

v[t-{2i  +  l)T,kTs]  =  0,for  [(2i-h2)iV  +  A:]T^  <0 
{0  <  t<  2NTs)  ,  (t  <  -2)  . 


Thus,  the  set  of  significant  i  =  (0,  — 1). 
Case  2i  k>  N 
For  this  set  of  k  values, 


p[t-2iT,kTs\  = 

(0  < 


0  ,  for  2iNT^  >  2N% 
t  <  2NT,)  ,  (i  >  1) 


and 


p[t-  {2i  + 1)7’,  kTs]  =  0  ,  for  [(2i  +  2)N  +  fc]  T,  <  0 

+ 

int 

Thus,  the  set  of  significant  i  =  ^0,  —  1,  •  •  • ,  —  • 

Case  3:  —N  <  A:  <  0 
For  this  set  of  k  values, 

p[t-  2iT,  kTs]  =  0  ,  for  {2iN  +  k)Ts  >  2NTs 
{0  <  t<  2NTs)  ,{i>2) 

and 


(0  <  t  <  2NTs)  ,  i<- 


k 

2N 


p[t-(2i+l)7’,A;Ts]  =  0  ,  for  (2i  +  2)JVrs  <  0 
(0  <  t  <  2NTs)  ,  (i  <  -1)  . 

Thus,  the  set  of  significant  i  =  (0, 1). 

Case  4z  k  <  —N 
For  this  set  of  k  values. 


and 


p[t-  2iT,  A:T*]  = 

(0  < 


0  ,  for  {2iN  +  k)Ts  >  2NT 
k 

'2N 


t  <  2NTs) 


-  I 


1  + 
int 


p[t-  {2i  +  1)7’, kTs]  =  0  ,  for  {2i  +  2)NTs  <  0 
(0  <  t  <  2NTs)  ,  (i  <  -1)  . 

Thus,  the  set  of  significant  i=  ^0, 1,  •  -  - ,  [~^]^j)- 

APPENDIX  D;  The  Four  Integrals  of  Equation  (28) 

Case  1:  0  <  k2  <  ki  —  N  —  1 


(■^1:^1  + 
L  2iV 


^  -2T  L  21V 

Fi(fci,A:2)  =  — jf  n  (•) 


dt 


(Dl) 
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An  exact  evaluation  of  ^2)  seems  to  be  unreasonably  complicated  to  carry  out.  Prom  the  properties  of  p  [t,  hTs\ 

(Appendix  B),  the  following  approximation  is  made. 


p[t-2iT,kT,]  “  p[t-{2i  +  l)T,kT,]^-^ 


Each  factor  in  (•)  (see  (29))  then  becomes 


(for  0  <  t  <  2NTs)  . 


1  sin  47r/io  1  sin  47r/ii 


4  sinTr/io  4  sinTr/ij 
=  cos  7r ho  cos  irhi  cos  27r/io  cos  27r/ii  . 

The  last  step  follows  from  (A6).  Fi(A;i,fc2)  then  becomes 

(fci ,  Ala)  —  [cos irho  cos  Trhi  cos  2itho  cos  2‘Khi^'^  [  int  . 

Case  2:  fci  —  AT  <  A12  <  fci  -  1 


(D2) 


(D3) 


(D4) 


1  ^ 

F2(Ai,A|2)  =  — jf  ][J(.)tft 


For  this  case,  the  integrand  can  be  shown  to  be 

1 
2 


^  (cos 3/3i  +  cos  /?i)  •  ^  (cos  3/^2  +  cos  ^2) 


+1  (cos  3/33  +  cos  ^^3)  •  i  (cos  3/?4  +  cos  /?4)  , 


using  (A6),  where 


/?!  =  2-Khop  [t,  ATs] 

/?2  =  27rhip  [t  -  T,  kTs] 

^3  =  2Trhop[t-2T,kT,] 

/34  =  2Trhip[t-3T,kT,]  . 

These  are  shown  below,  using  Appendix  B,  for  the  range  of  the  integral  {0,2NTs): 


A. 


A 

2xt^ 


A 


A 

2a’/), 


(D5) 


(D6) 
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Thus,  the  integral  becomes  the  sum  of  four  parts: 


1 


2NTs 


f  fNT, 

[I  * 


rNT,  r(2N-^k)T^  p2NT^  ^ 

J(N+k)n  JnT^  J{2N-irk)%  J 


/(N+k)T^  JNT, 

Evaluating  these  four  and  summing  produces,  after  many  tedious  steps,  the  result  is 

{^2  -  fei) 


Sirho  (k2  —  ki)  Tvho  (k2  ~  ki) 

- - - +COS - - 


0/356  3i  k\  1  ^  ^  TV  ~|“  k\ 


Zirhi  (k2  —  fci )  irhi  (k^  —  ki) 

4-  cos - - -  +  cos - - - 

N  N 


j_  rr _ ^ 

■47r  U3(/i§- 


+  ■ 


3hi 


+ 

+ 

+ 


[9/ig  -hi  hi-  hi 

ho  .  ho 


hi 


hi)  h^-9hl\ 

hi 


,  Zirhi  (fe  —  ki) 
N 


+ 


[9ftf  -hi  hi -hi 
ho  Zho 


.  irhi  (kz  —  fci) 

N 

.  irho  (kz  -  ki) 

TV 

.  Zithoikz-ki) 
- - N - 


}■ 


1  ® 

Hkukz)  =  ^]^ 

We  follow  the  same  procedure  used  to  obtain  Fz{ki,kz)',  accept  that  now  we  need 

/?!  =  2-Khop  [t,  fci;] 

02  =  2-Khip  [t  —  T,  AiT,] 

/?3  =  27rftoP  [t  +  2T,  fer,] 

/?4  =  27r/iip[t +  r,fcT,]  . 

These  are  shown  below  for  the  range  (0,2iVT,): 


A. 

^0 


NT,  (N  +  kyr,  2NT, 


A 

2;r^ 


NT,  (N+kyr,  2NT, 


(D7) 


(D8) 

(D9) 


A 

TxfiQ 


0^  2]^, 
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The  integral  becomes  the  sum  of  four  parts, 


2NT, 


requiring  many  tedious  steps  to  produce  the  result: 


I  {  fkT,  ^NT,  i‘(.N+k)T,  ^2NT,  \ 

JkT,  JnT.  J(N+k)T,) 


F3{kr,k2)  =  i 


J  _  (^2  -  fei)'  Sirho  {k2  -  ki)  ^  ^  tt/iq  (fc2  ~  ^i) 


N 


N 


N 


Sirhi  {k2  —  ki)  tt/ii  (ko  —  ki) 

+  - N - - N - 


+ 

+ 

OdS6  4j  iV  -|-  k\  1  ^  A^2  ^  -^-^0  —  1 


±/[ _ h. 

47r  \[3(hf- 


+ 


3fti 


/ig)^9h?-hgj 


.  Bit  hi  (k2  —  ki) 
- 


hi 

hi  1 

irhi  (^2  ■' 

-ki) 

M 

-9hl^ 

hf-hgj 

N 

ho 

ho 

Trho  (*2  - 

-ki) 

hi 

-Wi  ^ 

hg-h?J 

N 

[_ 

ho 

Bho 

.  37r/io 

{ka- 

[3  (/»§-*?) 

ml- 

hi. 

sin 

N 

1  ” 

n  (■) 

For  this  range  of  ^2  we  make  the  approximations, 


dt 


p[t-2ir,kT,]  “  p[t-(2i  +  l)T,Ar,]S- 
(forO  <  t<2NTs) 

and  follow  the  steps  used  in  obtaining  Fi  (ki,k2)  to  produce  the  result, 

-  ,  r fcg-fci  I  + 

F4(A:i,fc2)  —  [cosirhocos5rhicos27r/iocos27r/ii]  J<nt 


(DIO) 


(Dll) 


(D12) 


(D13) 
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